Inhomogeneous superconducting state in quasi-one-dimensional systems 
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We report on results of theoretical study of non-uniform superconducting states in quasi-one- 
dimensional systems, with attractive interactions and Zeeman splitting between electron spins. Using 
bosonization to treat intrachain electron-electron interactions, and a combination of renormalization 
group and mean-field approximation to tackle interchain couplings, we obtain the phase diagram 
of the system, and show that the transition between the uniform and non-uniform superconducting 
phases is a continuous transition of the commensurate-incommensurate type. 
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The possibility of a superconducting state with inho- 
mogeneous order parameter, stabilized by a sufficiently 
large Zeeman splitting between electrons with oppo- 
site spin orientations due to either an external mag- 
netic or internal exchange field, was suggested more 
than thirty years rfgo by Fulde and Ferrelfl and Larkin 
and Ovchinnikov.B Since then this Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state has been the subject of a num- 
ber of theoretical studies, but no direct evidence of its 
existence has ever been found in conventional supercon- 
ductors. More recently it has attracted renewed inter- 
est in the context of organic* pbeavy-fermion, and high- 
T c cuprate superconductors, Bc3 as these new classes of 
superconductors are believed to provide conditions that 
are favorable to the formation of FFLO state due to their 
quasi-one or two-dimensionality as well as unconventional 
pairing symmetry. Indeed, some gynerimjental evidence 
of its existence has been reported.BE3li3Ej 

The following picture emerged from the early theoreti- 
cal studies (mostly of mean-field type) of conventional 
s-wave superconductors subject to a Zeeman field B. 
For sufficiently high field the system is in the normal 
state. As the field strength decreases, at low tempera- 
tures the system undergoes a second-order phase tran- 
sition at B = B C 2(T) into the FFLO superconducting 
state. As the field strength further decreases, another 
phase boundary is encountered at B c i{T), and the sys- 
tem goes through another phase transition into the usual 
BCS superconducting state with uniform superconduct- 
ing order parameter. While it is much more difficult to 
locate the position of B c \ (T) than B C 2 (T) (even in mean- 
field theory), as well as to address the nature of the tran- 
sition there, it has been widely assumedc3 that this is a 
first-order phase boundary, across which the momentum 
of the order parameter and the magnetization change dis- 
continuously. This viewpoint was disputed in Ref. ^, in 
which the authors argue that the transition at B C \(T) is 
of second order. Thus the nature of this transition is an 
unsettled issue. 

In this paper we study quasi-one-dimensional (Q1D) 
superconductors subject to a Zeeman field, and the possi- 
bility of formation of FFLO states in these systems. Our 
motivation comes from two considerations. First of all, 
some of the experimental candidates for FFLO state are 



made of weakly coupled chains and therefore Q1D. Sec- 
ondly, it is known in that fluctuations are much stronger 
in low-dimensional systems than in 3D systems, and 
mean-field theories are much less reliable there.cJ On the 
other hand the non-perturbative machinery developed 
for studying one-dimensional interacting electron systems 
(especially bosonizationu) allows us to go beyond mean- 
field theory and treat the intrachain electron-electron 
correlation exactly in Q1D systems. In this paper we 
will take an approach that is similar to the one used in 
Ref. |2(| namely to treat the intrachain electron-electron 
interaction exactly using bosonization, and tackle the in- 
terchain couplings using a combination of renomalization 
group (RG) analysis and mean-field approximation. Us- 
ing this approach we are able to make a number of quan- 
titative and reliable predictions about the FFLO state 
in these systems. In particular, we will show that the 
phase transition at B c \ is continuous in these systems, 
and work out its critical properties. For the sake of sim- 
plicity and concreteness, we restrict our discussion to zero 
temperature throughout the paper. 

We start by considering a one-dimensional electron gas 
with attractive interactions. In the bosonized form, the 
Hamiltonian reads 



H — H c + H s + H z 



(1) 



where H c and H s are the Hamiltonian for the charge 
and spin sectors (which are decoupled, signaling the spin- 
charge separation) :l3 



H a = I dx\-^ 



+ V a COS(V87T0 Q ) f , 

(2) 



where a = c or s, and H z is the Zeeman coupling: 



H z = g^ B BSl ot = y -^9^bB J dxd x (f> s (x) 



(3) 



In these equations C and S are bosonic charge and spin 
fields related to the (coarse-grained) charge and spin den- 
sities: 

[2 FT 
Pi x ) = \ -d x (f) c (x), S z (x) = \ — d x <f> s (x); (4) 
V 7T V 2ir 
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while 9 a are their dual fields satisfying 

[(/) a (x),d x >0 a (x')} = iS(x-x'). 



B = B cl = 2 A s / 'guB, 



(10) 



(5) 



For attractive interactions, we typically have the Lut- 
tinger liquid parameters K c > 1 and K s < 1 (for non- 
interacting electrons, we have K c = K s = 1). If the 
ID electron gas is sufficiently far away from lattice com- 
mensuration, which we assume to be the case here, V c 
(which measures the strength of 4k f Umklapp scatter- 
ing) may be set to zero. Thus H c takes the form of free 
massless bosons. On the other hand, in the spin sector 
H s has the form of 1+1D quantum sine-don model, and 
for K s < 1, V s (which measures the strength of back 
scattering between electrons with opposite spins) is rele- 
vant in the RG sense; at low energies it opens up a gap 
A s ~ v„A\y,/v B A 3 ] 1 ^- aK ^ (A is the ultraviolet cutoff) 
for spin excitations .c3 The elementary spin excitations 
are massive solitons (kinks and anti- kinks) of the cf> s field, 
which carry spin ±1/2H3 This spin gap A s is the ana- 
log of quasi-particle gap in the BCS theory of higher di- 
mensional superconductors. The fundamental difference 
here, however, is that in ID there is no long-range su- 
perconducting order; instead the correlation function of 
the Cooper pair operator decays with a power law. The 
power law exponent can be calculated using the explicit 
representation of electron operators in terms of boson 
fields: 



i>\,<, = N a exp[i\kfx - i^x^ia 



(6) 



where A = ±1 represents left/right movers, a — ±1 rep- 
resents up/down spin particles, N a is the Klein factor 
that also includes a normalization constant, and 

$ A , CT = VV2[(#c - A0 C ) + a{9 s - \cp s )]. (7) 
Thus the singlet pair correlation function (at T = 0) 

(V'+i + i(z)V'ii-i(z)V'-i-i(2 y )V'+i+i(20) « 

(exp[iV2^(9 c (x) - c (a/))]}(exp[iV^Ffa.(aO - ^{x'))]) 

(8) 
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where the scaling dimension 



1 

2RZ 



(9) 



Here we have used the fact that the spin field 4> s (x) is 
long-range ordered in the spin-gapped phase. 

Let us now consider the effect of Hz- In Hz 
the Zeeman field couples to the soliton density and 
plays the role of chemical potential of spin solitons. 
In fact, H s + Hz takes-,,&xactly the form of the 
Pokrovsky-Talapov modelEafij which was introduced 
to study the two-dimensional classical commensurate- 
incommensurate (CIC) transition. In our context, we 
thus expect a continuous CIC transition at 



beyond which spin solitons start to proliferate in the 
ground state. Eq. |l^ is an exact result because the Zee-, 
man field couples to S 1 * * which is a conserved quantity.EJ 
In the incommensurate phase, the spin solitons form a 
spinless Luttinger liquid with its own bosonized Hamil- 
tonian, which describes the low-energy spin excitations 
of the system: 

In the long-wave length limit, the soliton density field 
(frsoiix) is related to the spin field <p s through 



S (x) = cj) so i(x)/V2- 



l{B)x + const., (12) 



where n so i is the soliton density of the ground state. In 
the limit B — > B c + + , the solitons become extremely 
dilute and the repulsive interaction among them become 
irrelevant; they can be treated as spinless free fermions.E2l 
As a consequence of this we have (i) K so i = 1 and (ii) 
n S oi(B) cc (B — B c ) x / 2 in this limit. Using these results 
we find in the incommensurate phase the superconduct- 
ing correlation function 



(ip f +1+1 {xy J l 1 _ 1 (x)^- 1 - 1 (x')^ +1+1 {x')) 
tx exp[iQ(B)(x - x')]\x - x'\~ 2 ^- 



(13) 



where Q(B) = im so i{B); approaching the phase bound- 
ary: B -> B c + 0+, we have Q{B) oc (B - B c )^ 2 and 



w 2K r 4 



= 6 



(14) 



This incommensurate phase (in the spin sector) is the 
ID analog of the FFLO phase in higher dimensional sys- 
tems, as the appearance of the spin solitons in the ground 
state induces an oscillatory phase in the superconduct- 
ing correlation function, Eq. (13). Also the ground state 
now has a finite magnetization as in the FFLO phase, 
and the additional fluctuation due to the soliton liquid 
makes the superconducting correlation function decays 
faster in the incommensurate phase, in (loose) analogy 
to the fact that appearance of unpaired quasiparticles 
reduces the size of the superconducting order parameter 
in the FFLO phase. We emphasize again that here there 
is no long-range superconducting order in either the com- 
mensurate or incommensurate phases; also the CIC tran- 
sition is continuous as n so i increases continuously from 
zero as B crosses B c i; both the magnetization and wave 
vector of oscillation Q are proportional to n so i. 

We now turn to the discussion of interchain couplings. 
The three leading potentially relevant perturbations to 
the decoupled Luttinger liquid fixed (dLL) point are 
single electron hopping H e , Cooper pair hopping (or 
Josephson tunnelpj|g) Hj, and interchain 2k / back scat- 
terings H C / SDW £3 For attractive interactions (K c > 1), 
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Hc/sdw is less relevant than Hj, in both the commen- 
surate and incommensurate phases. H e is irrelevant in 
the commensurate phase due to the presence of a spin 
gap. Since in this case the scaling dimension for Hj is 

O = 2U = l/K c < 2, (15) 

we conclude Cooper pair hopping is the leading relevant 
perturbation at the dLL fixed point, and the system flows 
toward a superconducting phase with long-range super- 
conducting order once interchain coupling is turned on 
in the commensurate phase. 

Now let us consider the incommensurate phase. Right 
after the system enters the incommensurate phase (B — > 
B c i + + ), we have 

Z'j = 2? ac = l/K c +l/2<2, (16) 

thus Hj is still relevant, albeit having a higher scaling 
dimension than that in the commensurate phase. How- 
ever in this case H e may also be relevant, as there is no 
longer a spin gap in this case. We find in this case the 
scaling dimension of H e to be 

C = \{K c + l/K c ) + b -. (17) 

We thus find that for K c > 3/2, £' 7 < C e , and Hj is 
the leading relevant perturbation at the dLL fixed point 
which drives the system to the Q1D superconducting 
FFLO phase once interchain coupling is turned on. On 
the other hand, for 1 < K c < 3/2, £' e < £' 7 < 2, and H e is 
the leading relevant perturbation at the dLL fixed point; 
in this case the system flaws toward the high-dimensional 
Fermi- liquid fixed point B3 These results are summarized 
in a schematic phase diagram, Fig. [l| The phase bound- 
ary separating the Fermi liquid and the two supercon- 
ducting phases are likely to be first-order since they are 
determined by the crossing of the scaling dimensions of 
two different relevant operators at the dLL fixed point; 
on the other hand as we will argue below, the transition 
from uniform to FFLO superconducting phases is contin- 
uous. We emphasize in this phase diagram we assume the 
Zeeman field B is not too strong; if the Zeeman splitting 
is so strong as to be comparable to, say the Fermi en- 
ergy, the continuum Luttinger liquid description of Q1D 
systems breaks down. 

To address the nature of the transition between uni- 
form and FFLO superconducting phase, we focus on 
the pair hopping process and neglect other perturbations 
that are less relevant: 

Hj=-ijY^ I rf^I^+i+i^-i-i^-i-iV'+i+i 

(ij) 

+ V'Vi-i^-i+iV'-i+iV'+i-i + h.c] 

= ~tj E / dxcos[V2^(9l - e{)] cos[VW s - (18) 

(ij) 



where i and j are chain indices, tj is the pair hopping 
matrix element (or Josephson coupling strength), (ij) 
stands for neighboring chains, tj oc tj, and h.c. stands 
for Hermitian conjugate. In the case of decoupled Lut- 
tinger liquids, there is spin-charge separation and the 
CIC transition occurs in the spin sector. As we see in Eq. 
( |l8| ) , interchain pair hopping couples the spin and charge 
fields. On the other hand since the system is in the su- 
perconducting phase (uniform or non-uniform) in which 
the charge field 8 C is long-range ordered, in studying the 
transition driven by B we may use a mean-field approx- 
imation and replace cos[v / 27r(^c — in Eq. ( |l8| ) by 
its expectation value: (cos[V2tt(9' 1 c — 8 J C )]) = C. Clearly 
this expectation value depends on B and it will also de- 
velop a dependence on x in the incommensurate phase; 
however as long as the dependence is smooth across the 
transition (which would be the case if the transition is 
continuous as we will show to be the case), we can treat 
it as a constant. Thus in the mean-field approximation 
Hj becomes 

Hf F = -Ctj [ Gfecos[V2^ s - 4)]. (19) 

(ij) 

Eq. ( jl9| ) can also be obtained more formally by integrat- 
ing out the fluctuations of the 6 C field on top of its ex- 
pectation value in the Lagrangian formalism, which will 
yield a slightly renormalized coupling C. The quantum 
Hamiltonian of H s + Hz + Hj IF can be mapped onto 
the problem of classical CIC transition driven by B at 
finite temperatures, in d + 1 dimensions (d is the phys- 
ical dimension of the quantum problem we study here). 
It is known that the CIC transition in higher dimensions 
is still continuous, but the critical behavior is very dif- 
ferent from the d = 1 case considered earlier; in this case 
the density of domain walls (that consist of solitons of 
individual chains aligned with true long-range order) de- 
pends logarithmically on the distance from criticalityS 

n wall ^l/\og{\B-B c \- 1 ) (20) 

as B — > B c + + . The wave vector of the inhomogeneous 
superconducting order parameter Q and the magnetiza- 
tion are both proportional to n wa u and thus have the 
same dependence on B near criticality. We note that 
while we obtained these results by making a mean-field 
approximation to the (long-range ordered) charge fields, 
the main conclusion that the transition is continuous 
should be robust; this follows simply from the fact that 
the domain walls (whose appearance drives the transi- 
tion) repel each other, which is clearly the case here. 
The logarithmic dependence of n wa u on B — B c then fol- 
lows from the exponentially weak repulsion between the 
domain walls. These in turn justify the validity of the 
mean-field approximation employed. 

To summarize, we studied formation of non-uniform 
superconducting state in quasi-one-dimcnsional systems. 
Among our results include a phase diagram in terms of 
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the Zeeman field and Luttinger liquid parameter. We 
also showed that the transition between the uniform and 
non-uniform superconducting states is continuous. 

The author has benefited greatly from discussions with 
David Huse. This work was supported by NSF DMR- 
9971541 and the A. P. Sloan Foundation. 



1 P. Fulde and A. Ferrell, Phys. Rev. 135, A550 (1964). 

2 A. I. Larkin and Yu. N. Ovchinnikov, Sov. Phys. JETP 20, 
762 (1965). 

3 K. Gloos et al, Phys. Rev. Lett. 70, 501 (1993). 

4 G. Yin and K. Maki, Phys. Rev. B 48, 650 (1993). 

5 M. R. Norman, Phys. Rev. Lett. 71, 3391 (1993). 

6 H. Burkhardt and D. Rainer, Ann. Physik 3, 181 (1994). 

7 H. Shimahara, Phys. Rev. B 50, 12760 (1994). 

8 G. Murthy and R. Shankar, J. Phys. Condens. Matter 7, 
9155 (1995). 

9 N. Dupuis, Phys. Rev. B 51, 9074 (1995). 

10 R. Modler et al, Phys. Rev. Lett. 76, 1292 (1996). 

11 M. Tachiki et al, Z. Phys. B 100, 369 (1996). 

12 P. Gegebwart et al, Ann. Physik 5, 307 (1996). 

13 H. Shimahara, S. Matsuo, and K. Nagai, Phys. Rev. B 53, 
12284 (1996). 

14 K. Maki and H. Won, Czech. J. Phys. 46, 1035 (1996). 

15 K. V. Samohkin, Physica C 274, 156 (1997). 

16 A. I. Buzdin and H. Kachkachi, Phys. Lett. A 225, 341 
(1997). 

17 K. Yang and S. L. Sondhi, Phys. Rev. B 57, 8566 (1998). 

18 W. E. Pickett, R. Weht, and A. B. Shick, Phys. Rev. Lett. 
83, 3713 (1999). 

19 K. Yang and D.F. Agterberg, Phys. Rev. Lett. 84, 4970 
(2000); D. F. Agterberg and K. Yang, |cond-mat/0006344 . 

20 S. Manalo and U. Klein, J. Phys. Condens. Matter 12, L471 
(2000). 

21 J. Singleton, J. A. Symington, M.-S. Nam, A. Ardavan, M. 
Kurmoo, and P. Days, J. Phys. Condens. Matter 12, L641 
(2000). 

22 D. F. Agterberg and K. Yang, |cond-mat/0006344 

23 See, e.g., D. Saint- James, G. Sarma and E. J. Thomas 
(eds.), Type II Superconductors, Pergamon Press, New 
York (1969); L. W. Gruenberg and L. Gunther, Phys. Rev. 
Lett. 16, 996 (1966); as well as Refs. (?]-§[ 

24 As a matter of fact in the past a number of authors have 
studied FFLO state in superconductors in the ID limit us- 
ing mean-field theory and a strictly ID electron dispersion; 
see A. I. Buzdin and V. V. Tugushev, Sov. Phys. JETP 
58, 428 (1983); K. Machida and H. Nakanishi, Phys. Rev. 
B 30, 122 (1984); A. I. Buzdin and S. V. Polonskii, Sov. 
Phys. JETP 66, 422 (1987). In these mean-field theories 
one assumes a long-range ordered superconducting order 
parameter. However it is known (also see later) that quan- 
tum as well as thermal fluctuations do not allow long-range 
superconducting order in ID, at zero or finite temperature, 
and thus put the validity of mean-field theory in question 
in ID. 



A. O. Gogolin, A. A. Nersesyan and A. M. Tsvelik, 
Bosonization and Strongly Correlated Systems, Cambridge 
University Press, Cambridge (1998). 

' E. W. Carlson, D. Orgad, S. A. Kivelson, and V. J. Emery, 
Phys. Rev. B 62, 3422 (2000). 

For K s < 1/2, there also exist more complicated excitations 
that are bound states of kinks and anti-kinks. However it is 
unlikely that short-range attractive interactions can place 
the system in this regime (for a recent discussion of this 
point, as well as what may happen w hen K 3 becomes ver y 



small, see E. Papa and A. M. Tsvelik, cond-mat/0006422). 
Also as we will see later, the Zeeman field will only generate 
one type of solitons in the ground state. We thus do not 
consider complications in this regime here. 
V. L. Pokrovsky and A. L. Talapov, Phys. Rev. Lett. 42, 
65 (1979). 

H. J. Schulz, Phys. Rev. B 22, 5274 (1980). 

S. Sachdev, T. Senthil and R. Shankar, Phys. Rev. B 50, 

258 (1994). 

It was emphasized recently (V. J. Emery, E. Fradkin, S. 
A. Kivelson and T. C. Lubensky, Phys. Rev. Let t. 85 , 
2160 (2000); A. Vishwanath and D. Carpenter, cond 
mat/0003036) that interchain forward scattering are ex 



actly marginal operators at the dLL fixed point and can 
have important effects on the scaling dimensions of other 
perturbations. We do not consider these issues here. 
In principle there may be residual attractive interaction 
among quasiparticles that can destabilize the Fermi liquid 
fixed point that the system flows toward, and trigger a su- 
perconducting instability. An analysis of such possibilities 
is beyond the scope of the approach used here since we 
can only analyze the possible instabilities of the dLL fixed 
point in this approach. 

M. E. Fisher and D. S. Fisher, Phys. Rev. B 25, 3192 
(1982). 



B 





Nonunifrom (FFLO) 




Superconductor / 


Fermi Liquid 


/ 


Uniform Superconductor 



1.0 



1.5 

K 



2.0 



FIG. 1. Schematic phase diagram of coupled Luttinger 
liquids subject to a Zeeman field. The solid lines are first 
order phase boundaries while the dashed line is a second-order 
phase boundary. 
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